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Abstract. We study of the effect of turbulence on diffusion processes within magnetized medium. 
While we exemplify our treatment with heat transfer processes, our results are quite general and are 
applicable to different processes, e.g. diffusion of heavy elements. Our treatment is also applicable 
to describing the diffusion of cosmic rays arising from magnetic field wandering. In particular, we 
find that when the energy injection velocity is smaller than the Alfven speed the heat transfer is 
partially suppressed, while in the opposite regime the effects of turbulence depend on the intensity 
of driving. In fact, the scale Ia at which the turbulent velocity is equal the Alfven velocity is a new 
important parameter When the electron mean free path X is larger than Ia, the stronger the the 
turbulence, the lower thermal conductivity by electrons is. The turbulent motions, however, induces 
their own advective transport, that can provide effective diffusivity. For clusters of galaxies, we find 
that the turbulence is the most important agent for heat transfer. We also show that the domain of 
applicability of the subdiffusion concept is rather limited. 
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ASTROPHYSICAL TURBULENCE 

The magnetohydrodynamic (MHD) approximation is widely used to describe the actual 
magnetized plasma turbulence over scales that are much larger than both the mean 
free path of the particles and their Larmor radius (see Kulsrud 2004). The theory 
of MHD turbulence has become testable recently due to numerical simulations (see 
Biskamp 2003) which confirm (see Cho & Lazarian 2005 and ref. therein) the prediction 
of magnetized Alfvenic eddies being elongated along magnetic field (see Shebalin, 
Matthaeus & Montgomery 1983, Higdon 1984) and provided results consistent with the 
quantitative relations for the degree of eddy elongation obtained in Goldreich & Sridhar 
(1995, henceforth GS95). 

GS95 model assumes the isotropic injection of energy at scale L and the injection 
velocity equal to the Alfven velocity in the fluid Va, i.e. the Alfven Mach number 
Ma = {5V /Va) = 1. This model can be easily generalized for both Ma > 1 and Ma < 1 
at the injection. Indeed, if Ma > 1, instead of the driving scale L for one can use another 
scale, namely Ia, which is the scale at which the turbulent velocity gets equal to Va- For 
Ma ^ 1 magnetic fields are not dynamically important at large scales and the turbulence 
follows the incompressible Kolmogorov cascade V/ ~ /^/^ over the range of scales [L, Ia] ■ 
This provides I a ~ L^M^^ . If Ma < 1, the turbulence obeys GS95 scaling (also called 
"strong" MHD turbulence) not from the scale L, but from a smaller scale krans ~ 
(Lazarian & Vishniac 1999), while in the range [L, hrans] the turbulence is "weak". 



BASICS OF DIFFUSION IN MAGNETIZED PLASMAS 



The issue of diffusion in magnetized plasma has been mostly dealt with in the context of 
heat transfer. Heat transfer in turbulent magnetized plasma is an important astrophysical 
problem which is relevant to the wide variety of circumstancies from mixing layers in 
the Local Bubble (see Smith & Cox 2001) and Milky way (Begelman & Fabian 1990) 
to cooling flows in intracluster medium (ICM) (Fabian 1994). The latter problem has 
been subjected to particular scrutiny as observations do not support the evidence for 
the cool gas (see Fabian et al. 2001). This is suggestive of the existence of heating that 
replenishes the energy lost via X-ray emission. Heat transfer from hot outer regions is 
an important process to consider in this context. 

It is well known that magnetic fields can suppress thermal conduction perpendicular 
to their direction. The issue of heat transfer in realistic turbulent magnetic fields has 
been long debated (see Bakunin 2005 and references therein). An influencial paper by 
Narayan & Medvedev (2001, henceforth NMOl) obtained estimates of thermal conduc- 
tivity by electrons using the GS95 model of MHD turbulence with the velocity Vi at the 
energy injection scale L that is equal to the Alfven velocity Va, i.e. the turbulence with 
the Alfven Mach number Ma = (Vl/Va) = 1. This is rather restrictive, as in the ICM 
Ma > 1, while in other astrophysical situations Ma < 1. Below we discuss turbulence 
for both Ma > 1 and Ma < 1 and compare the particle diffusion to that by turbulent fluid 
motions. For heat transfer, we refer to our results in Lazarian (2006, 2007). 

Let us initially disregard the dynamics of fluid motions on diffusion, i.e. consider dif- 
fusion induced by particles moving along static magnetic fields. Magnetized turbulence 
in the GS95 model is anisotropic with eddies elongated along (henceforth denoted by ||) 
the direction of local magnetic field. Consider isotropic injection of energy at the outer 
scale L and dissipation at the scale l±,mim where ± denotes the direction of perpendicular 
to the local magnetic field. NMOl observed that the separations of magnetic field lines 
for ro < l±^min mostly influenced by the motions at the scale l±,mim which results 
in Lyapunov-type growth: ~ roexp(Z//|| „„■„). This growth is similar to that obtained in 
earlier models with a single scale of turbulent motions (Rechester & Rosenbluth 1978, 
Chandran & Cowley 1998). This is not surprising as the largest shear that causes field 
line divergence is provided by the marginally damped motions at the scale around l±^min- 
In NMOl ro is associated with the size of the cloud of electrons of the electron Larmor 
radius riar,particie- They find that the electrons should travel over the distance 

l-RR ~ ^||,m/;3ln(Z^^„„>,/ riar,particle) (1) 

to get separated by 

Within the single-scale model which formally corresponds to L = /y = l±,min 
the scale Lrr is called Rechester-Rosenbluth distance. For the ICM parameters the 
logarithmic factor in Eq. ([T]) is of the order of 30, and this causes 30 times decrease 
of thermal conductivity for the single-scale model^. In the multi-scale models with a 
limited (e.g. a few decades) inertial range the logarithmic factor stays of the same order 



^ For the single-scale model Lrr ^ 30L and the diffusion over distance A takes Lrr/L steps, i.e. ^ 
LrrL, which decreases the corresponding diffusion coefficient K'partic/e.img/f ^ 

A^/5t by the factor of 30. 



but it does not affect the thermal conductivity, provided that L ^ / 1 Indeed, for the 
electrons to diffuse isotropically they should spread from riar.partide to ^- The GS95 
model of turbulence operates with field lines that are sufficiently stiff, i.e. the deviation 
of the field lines from their original direction is of the order unity at scale L and less for 
smaller scales. Therefore to get separated from the initial distance of to a distance 
L (see Eq. ([5]) with Ma = 1), at which the motions get uncorrelated the electron should 
diffuse the distance slightly larger (as field lines are not straight) than VlL, which is 
much larger than the extra travel distance ~ 30/|| Explicit calculations in NMOl 
support this intuitive picture. 

DIFFUSION FOR Ma > 1 

Turbulence with Ma > 1 evolves along hydrodynamic isotropic Kolmogorov cascade, 
i.e. Vi ~ Vl{1/L) over the range of scales [L, Ia], where 

lA-L{VA/VL)^=LMf, (2) 

is the scale at which the magnetic field gets dynamically important, i.e. V/ = Va. This 
scale plays the role of the injection scale for the GS95 turbulence, i.e. V/ ~ 1^ (^±/^a) ^'^^^ 
with eddies at scales less than Ia geting elongated in the direction of the local magnetic 
field. The corresponding anisotropy can be characterized by the relation between the 
semi-major axes of the eddies 

l\\^L{ljLfl^M^\ Ma>1. (3) 

where || and ± are related to the direction of the local magnetic field. In other words, for 
Ma > 1, the turbulence is still isotropic at the scales larger to Ia, but develops (/^/ Ia) 
anisotropy for I < Ia- 

If particles (e.g. electrons) mean free path A ^ Ia, they stream freely over the distance 
of I A- For particles initially at distance l±^min to get separated by L the required travel 
is the random walk with the step Ia, i.e. the mean-squared displacement of a particle 
till it enters an independent large-scale eddy ~ 1\{L/Ia), where L/Ia is the number 
of steps. These steps require time 5t ~ {L/ 1 a)Ia/ C\v particle, where Vpartide is electron 
thermal velocity and the coefficient Ci = 1/3 accounts for ID character of motion along 
magnetic field lines. Thus the electron diffusion coefficient is 

^particle = / 5t ^ {\ / 3) I aV particle, Ia < ^ , (4) 

which for ^ A constitutes a substantial reduction of diffusivity compared to its 
unmagnetized value Kunmagn = ^Vpartide- We assumed in Eq. (HI) that L 30/|| (see 
§2.1). 

For X <^ Ia <^ L, Kpartide ~ ^/^iKunmagn as both the Lrr and the additional distance 
for electron to diffuse because of magnetic field being stiff at scales less than Ia are 
negligible compared to L. For Ia L, when magnetic field has rigidity up to the scale 
L, it gets around 1/5 of the value in unmagnetized medium, according to NMOl. 



Note, that even dynamically unimportant magnetic fields do influence heat conduc- 
tivity over short time intervals. For instance, over time interval less than l\/CiKunmagn 
the diffusion happens along stiff magnetic field lines and the difference between paral- 
lel and perpendicular diffusivities is largely. This allows the transient existence of sharp 
small-scale temperature gradients. 



DIFFUSION FOR Ma < 1 

It is intuitively clear that for Ma < 1 turbulence should be anisotropic from the injection 
scale L. In fact, at large scales the turbulence is expected to be weaAH (see Lazarian 
& Vishniac 1999, henceforth LV99). Weak turbulence is characterized by wavepackets 
that do not change their Z||, but develop structures perpendicular to magnetic field, i.e. 
decrease l± . This cannot proceed indefinitely, however. At some small scale the GS95 
condition of critical balance, i.e. 1\\/Va ~ l±/Vi, becomes satisfied. This perpendicular 
scale Itrans Can be obtained substituting the scaling of weak turbulence (see LV99) 
Vi ~ Vi{l^/Ly/'^ into the critical balance condition. This provides In-ans ~ ^^'^ '•^^ 
corresponding velocity Vtrans ~ VlMa- For scales less than krans the turbulence is strong 
and it follows the scalings of the GS95-type, i.e. V/ ~ VL{L/l^y^/^M\^^ and 

/|| ~L(/i/L)2/3M//^ Ma<1. (5) 

For Ma < 1, magnetic field wandering in the direction perpendicular to the mean 
magnetic field (along y-axis) can be described by d{y^)/dx ~ (/)//|| (LV99), whereQ /y 

is expressed by Eq. ([5]) and one can associate l± with 2{y^) 

(y^^^^-^^^^A. l±< Itrans (6) 

For weak turbulence d{y^)/dx ~ LMa (LV99) and thus 

(/)l/2~Ll/2^1/2Mi, h> Itrans. (7) 

Fig. 1 confirms the correctness of the above scaling numerically. 

Eq. © differs by the factor M^ from that in NMOl , which reflects the gradual suppres- 
sion of thermal conductivity perpendicular to the mean magnetic field as the magnetic 
field gets stronger. Physically this means that for Ma < 1 the magnetic field fluctuates 
around the well-defined mean direction. Therefore the diffusivity gets anisotropic with 



^ The relation between the mean squared displacements perpendicular to magnetic field (y^) and the 
displacements x along magnetic field for x < can be obtained through the diffusion equation approach 

in §2.3 and Eq. ©. This gives (/)i/2 _ -^^mI 

^ The terms "weak" and "strong" turbulence are accepted in the literature, but can be confusing. As we 
discuss later at smaller scales at which the turbulent velocities decrease the turbulence becomes strong. 
The formal theory of weak turbulence is given in Galtier et al. (2000). 

^ The fact that one gets Z|| „„■„ in Eq. ^ is related to the presence of this scale in this diffusion equation. 
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FIGURE 1. Root mean square separation of field lines in a simulation of inviscid MHD turbulence, 
as a function of distance parallel to the mean magnetic field, for a range of initial separations. Each 
curve represents 1600 line pairs. The simulation has been filtered to remove pseudo-Alfven modes, which 
introduce noise into the diffusion calculation. From Lazarian, Vishniac & Cho 2004. 



the diffusion coefficient parallel to the mean field K\\^particie ~ l/3K'«nmag(i being larger 
than coefficient for diffusion perpendicular to magnetic field k± particle- 

Consider the coefficient K^^panwie for ^ 1- As NMOl showed, particles become 
uncorrected if they are displaced over the distance L in the direction perpendicular to 
magnetic field. To do this, a particle has first to travel Lrr (see Eq. ([1])), where Eq. dS]) 
relates /|| and l±jnin- Similar to the case in §2.1, for L 30/|| „,;,„ the additional travel 
arising from the logarithmic factor is negligible compared to the overall diffusion dis- 
tance L. At larger scales electron has to diffuse ~ L in the direction parallel to magnetic 
field to cover the distance of LM^ in the direction perpendicular to magnetic field di- 
rection. To diffuse over a distance R with random walk of LM^ one requires /L^M\ 
steps. The time of the individual step is L^/ K\^\^particie- Therefore the perpendicular diffu- 
sion coefficient is 

^L.partide = / {R^ /[KlparticleMl]) = KlpartideMl. Ma < 1 , (8) 

An essential assumption there is that the particles do not trace their way back over the 



individual steps along magnetic field lines, i.e. Lrr «L. Note, that for of the order 
of unity this is not accurate and one should account for the actual 3D displacement. This 
introduces the change by a factor of order unity (see above). 



TURBULENT DIFFUSIVITY 

Turbulent motions themselves can induce advective transport. In Cho et al. (2003) we 
dealt with the turbulence with Ma ~ 1 and estimated 

^dynamic ~ CdynLVt, > 1 , (9) 

where Cdyn ~ 0(1) is a constant, which for hydro turbulence is around 1/3 (Lesieur 
1990). If we deal with heat transport, for fully ionized non-degenerate plasma we assume 
Cdyn ~ 2/3 to account for the advective heat transport by both protons and electron^. 
Thus eq. ^ covers the cases of both Ma > 1 up to Ma ~ 1. For Ma < 1 one can 
estimate Kdynamic ~ (fio), where d is the random walk of the field line over the wave 
period ~ O)^ . As the weak turbulence at scale L evolves over time T ~ M^^(0^\ 
{y^) is the result of the random walk with a step d, i.e. (j^) ~ (T(o)(i^. According to 
eq.® and (|7]), the field line is displaced over time T by {y^) ~ LM^V^T. Combining 
the two one gets d^ ~ LMj^tO)^, which provides K^Jy^amic ~ Qyn^^L^l' which is 
similar to the diffusivity arising from strong turbulence at scales less than hrans, i-e- 
^dynamic ~ CdynhransVtrans- The total diffusivity is the sum of the two, i.e. for plasma 

Kdynam^c-{^/^)LVLMl Ma<\, (10) 

where /3 ~ 4. 



EXAMPLE: THERMAL CONDUCTIVITY 

In thermal plasma, electrons are mostly responsible for thermal conductivity. The 
schematic of the parameter space for Kpartide < ^dynamic is shown in Fig 2, where the 
the Mach number M^ and the Alfven Mach number Ma are the variables. For Ma < 1, 
the ratio of diffusivities arising from fluid and particle motions is Kdynamic/ '>^partide ~ 
PaMsMA{L/X) (see Eqs. ^ and (flOl)). the square root of the ratio of the electron to 
proton mass a = {me/nipy/^, which provides the separation line between the two re- 
gions m Fig. 2, /3aM, ~ (A/L)Ma. For \<Ma< [L/Xy/^ the mean free path is less 
than I A which results in Kpamde being some fraction of Kunmagn, while Kdynamic is given 
by Eq. Q. Thus Kdynamic / Kpartide ~ ^OiMs{L/X), i.e. the ratio does not depend on Ma 
(horisontal line in Fig. 2). When Ma > {L/Xy/^ the mean free path of electrons is con- 
strained by Ia- In this case Kdynamic / Kpartide ~ ^aMsM\ (see Eqs. dH) and @) . This 
results in the separation line [5aMs ~ M^^ in Fig. 2. 



This gets clear if one uses the heat flux equation q = —Kc\jT, where Kc = nkg Kiymimic/eiectr^ « is electron 
number density, and kg is the Bohzmann constant, for both electron and advective heat transport. 
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FIGURE 2. Parameter space for particle diffusion or turbulent diffusion to dominate: application to 
heat transfer. Sonic Mach number Ms is ploted against the Alfven Mach number Ma. The heat transport is 
dominated by the dynamics of turbulent eddies is above the curve and by thermal conductivity of electrons 
is below the curve. Here A is the mean free path of the electron, Lis the driving scale, and a = {me/mpYI'^, 
/3 w 4. The panel in the right upper coner of the figure illustrates heat transport for the parameters for a 
cool core Hydra cluster (point "F"), "V" corresponds to the illustrative model of a cluster core in EnssUn 
et al. (2005). From Lazarian (2007). 



HEAT TRANSFER IN INTRACLUSTER MEDIUM 



It is generally believed that Intracluster Medium (ICM) is turbulent. The considerations 
below can be used as guidance. In unmagnatized plasma with the ICM temperatures 
T ~ 10^ K and and density 10^-^ cm^-^ the kinematic viscosity r\unmagn ~ Vion^Um, where 
Vion and Xion are the velocity of an ion and its mean free path, respectively, would make 
the Reynolds number Re = LVi/rjunmagn of the order of 30. This is barely enough for 
the onset of turbulence. For the sake of simplicity we assume that ion mean free path 
coinsides with the proton mean free path and both scale as A ~ 3T^nZ\ kpc, where the 
temperature = kT /3 keV and n_3 = n/lO~^ cm~^. This provides A of the order of 
0.8-1 kpc for the ICM (see NMOl). 

It is accepted, however, that magnetic fields decrease the diffusivity. Somewhat 
naively assuming the maximal scattering rate of an ion, i.e. scattering every orbit (the 
so-called Bohm diffusion limit) one gets the viscosity perpendicular to magnetic field 
^7_L ~ VionrLar,ion, which is much Smaller than rfunmagn, provided that the ion Larmor 
radius rLar,ion <S hon- For the parameters of the ICM this allows essentially invicid mo- 



tion^ of magnetic lines parallel to each other, e.g. Alfven motions. 

In spite of the substantial progress in understading of the ICM (see EnBlin, Vogt 
& Pfrommer 2005, henceforth EVP05, EnBlin & Vogt 2006, henceforth EV06 and 
references therein), the basic parameters of ICM turbulence are known within the factor 
of 3 at best. For instance, the estimates of injection velocity Vl varies in the literature 
from 300 km/s to 10^ km/s, while the injection scale L varies from 20 kpc to 200 
kpc, depending whether the injection of energy by galaxy mergers or galaxy wakes is 
considered. EVP05 considers an illustrative model in which the magnetic field with the 
10 /iG fills 10% of the volume, while 90% of the volume is filled with the field of 5 ~ 1 
juG. Using the latter number and assuming Vl = 10^ km/s, L = 100 kpc, and the density 
of the hot ICM is 10^"^ cm^-^, one gets ^ 70 km/s, i.e. Ma > 1. Using the numbers 
above, one gets Ia ~ 30 pc for the 90% of the volume of the hot ICM, which is much less 
than Xion - The diffusivity of ICM plasma gets rj = viouIa which for the parameters above 
provides 7?e ~ 2 x 10-^, which is enough for driving superAlfvenic turbulence at the outer 
scale L. However, as Ia increases as 5^, Re gets around 50 for the field of 4 /iG, which 
is at the border line of exciting turbulence^. However, the regions with higher magnetic 
fields (e.g. 10 /iG) can support Alfvenic-type turbulence with the injection scale Ia and 
the injection velocities resulting from large-scale shear Vl{Ia/L) ~ VlM^^. 

For the regions of 5 ~ 1 /iG the value of I a is smaller than the mean free path of 
electrons A . According to Eq. (HI) the value of K^igctr is 100 times smaller than KspUzer- On 
the contrary, Kdynamic for the ICM parameters adopted will be ~ ^OKspUzen which makes 
the dynamic diffusivity the dominant process. This agrees well with the observations 
in Voigt & Fabian (2004). Fig. 2 shows the dominance of advective heat transfer for 
the parameters of the cool core of Hydra A ( 5 = 6 /iG, n = 0.056 cm L = 40 kpc, 
T = 2.7 keV according to EV06), point "F", and for the illustrative model in EVP05, 
point "V", for which B=l iiG. 

Note that our stationary model of MHD turbulence is not directly applicable to 
transient wakes behind galaxies. The ratio of the damping times of the hydro turbulence 
and the time of straightening of the magnetic field lines is ~ M^^. Thus, for Ma > 1, 
the magnetic field at scales larger than Ia will be straightening gradually after the hydro 
turbulence has faded away over time L/V^. The process can be characterized as injection 
of turbulence at velocity Va but at scales that increase linearly with time, i.e. as Ia + VAt. 
The study of heat transfer in transient turbulence and magnetic field "regularly" stretched 
by passing galaxies will be provided elsewhere. 



^ A regular magnetic field w {2mkTy l^c / [eX] that makes rLar,ion less than A and therefore rj^ < 
Vunmagn is just 10^^" G. Turbulent magnetic field with many reversals over ruir.itm does not interact 
efficiently with a proton, however As the result, the protons are not constrained until gets of the order 
of ruir.ioiT This happens when the turbulent magnetic field is of the order of 2 x 10^^(VL/10^km/s) G. At 
this point, the step for the random walk is ^ 2 x 10^^ pc and the Reynolds number is 5 x 10^". 
^ One can imagine dynamo action in which superAlfvenic turbulence generates magnetic field till Ia gets 
large enough to shut down the turbulence. 
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FIGURE 3. Energy density of compressive modes and Alfvenic slab-type waves, induced by CRs. The 
energy is transferred from the mean free path scale to the CR Larmor radius scale. If the mean free 
path falls below compressive motions cutoff or feedback suppression scale, the spectrum of slab waves 
becomes steeper. From Lazarian & Beresnyak (2006). 

ASPECTS OF COSMIC RAY DIFFUSION 

Diffusion of cosmic rays (CR) includes diffusion arising from CR scattering and that of 
wandering magnetic field lines. For our purposes, we disregard the diffusion perpendic- 
ular to magnetic field lines that arises from gyroresonance interactions. 

Naturally, the diffusion of magnetic field lines depends on the adopted model of 
turbulence. For instance, Lazarian & Beresnyak (2006) discuss a mechanism by which 
slab-type perturbations can be generated on the scale of CR gyroradius (seeFig. 3). These 
slab-type perturbations are likely to contribute subdominantly to field line wandering. 
Therefore below we concentrate, as in the rest of the paper, on the field wandering 
induced by Alfvenic modes of GS95 turbulence. 

Generalizing our arguments above in terms of the diffusion parallel and perpendicular 
to mean magnetic field, one can write that k^ ^t = '^||,ct(^^)^/(^z)^' where bR is a 
random walk step in the direction perpendicular to magnetic field. If the corresponding 
scale bz is less than Lrr, the motion along the field line is one dimentional diffusion that 
retraces its steps. In this case, 5z = {K.\^ „8t)^l'^ . If, following earlier authors (see Kota & 
Jokipii 2000, Web et al. 2006) we introduce a spatial field diffusion coefficient Kspat = 

{dR)'^/{dz), we easily get Ki KspatK}A^{dty/'^, which results in the subdiffusion in 



perpendicular direction, i.e. in the distance perpendicular to magnetic field growing as 
fc^cr^^ = KspaiK^^^^^{dty^^. Howcvcr, as soon as the distance to diffuse is much larger 
than Lrr (see Eq.©, the subdiffusive effects are negligible. 

Although CR velocities are of the order of light speed, for sufficiently small mean free 
paths the advection of CR by turbulent motions may become important. 

CONCLUDING REMARKS 

In the paper above we attempted to describe the diffusion by particle and turbulent 
motions for Ma < 1 and Ma > 1. Unlike earlier papers, we find that turbulence may 
both enhance diffusion and suppress it. For instance, when A gets larger than Ia the 
conductivity of the medium ~ M^^ and therefore the turbulence inhibits heat transfer, 
provided that Kpartide > ^dynamic- Aloug with the plasma effects that we mention below, 
this effect can, indeed, support sharp temperature gradients in hot plasmas with weak 
magnetic field. 

As discussed above, rarefied plasma, e.g. ICM plasma, has large viscosity for motions 
parallel to magnetic field and marginal viscosity for motions that induce perpendicular 
mixing. Thus fast dissipation of sound waves in the ICM does not contradict the medium 
being turbulent. The later may be important for the heating of central regions of clusters 
caused by the AGN feedback (see Churasov et al. 2001, Nusser, Silk & Babul 2006 
and more references in EV06). Note, that models that include both heat transfer from 
the outer hot regions and an additional heating from the AGN feedback look rather 
promissing (see Ruszkowkski & Begelman 2002, Piffaretti & Kaastra 2006). We predict 
that the viscosity for 1 jiG regions is less than for 10 jlG regions and therefore heating 
by sound waves (see Fabian et al. 2005) could be more efficient for the latter. Note, 
that the plasma instabilities in coUisionless magnetized ICM arising from compressive 
motions (see Schekochihin & Cowley 2006, Lazarian & Beresnyak 2006) can resonantly 
scatter particles and decrease A . This decreases further Kpanide compared to Kunmagn but 
increases Re. In addition, we disregarded mirror effects that can reflect electrons back 
(see Malyshkin & Kulsrud 2001 and references therein), which can further decrease 

^particle- 

All in all, we have shown that it is impossible to characterize the diffusion in mag- 
netized plasma by a single fraction of the diffusion coefficient of unmagnetized value. 
The diffusion depends on sonic and Alfven Mach numbers of turbulence and the corre- 
sponding diffusion coefficient may be much higher and much lower than the unmagne- 
tized one. As the result, turbulence can inhibit or enhance diffusivity depending on the 
plasma magnetization and turbulence driving. 

Our study indicates that in many cases related to ICM the advective heat transport by 
dynamic turbulent eddies dominates thermal conductivity. In addition, "subdiffusivity" 
is probably subdominant for many astrophysical problems. 

Acknowledgments The research is supported by the NSF Center for Magnetic Self 
Organization in Laboratory and Astrophysical Plasmas. Help by A. Beresnyak is ac- 
knowledged. 



REFERENCES 



1. Bakunin, O.G. 2005, Plasma Phys. and Contr. Fus. 47, 1857 

2. Begelman, M. & Fabian, A. 1990, MNRAS, 244, 26 

3. Biskamp, D. 2003, Magnetohydrodynamic Turbulence. (Cambridge: CUP) 

4. Chandran, B. & Cowley, S. 1998, Phys. Rev. Lett., 80, 3077 

5. Cho J., Lazarian A., 2005, Theoret. Comput. Fluid Dynamics, 19, 127 

6. Cho, J. & Lazarian, A. 2002, Phys. Rev. Lett., 88, 5001 

7. Cho, J., Lazarian, A., Honein, A., Kassions, S., & Moin, P 2003, ApJ, 589, L77 

8. Churazov, E., Bmggen, M., Kaiser, C, Bohringer, H., & Forman W. 2001, ApJ, 554, 261 

9. Galtier, S., Nazarenko, S., Newel, A. & Pouquet, A. 2000, J. Plasma Phys., 63, 447 

10. EnBlin, T., Vogt, C. & Pfrommer, C. 2005, in The Magnetized Plasma in Galaxy Evolution, Eds. K.T. 
Chyzy, K. Otminowska-Mazur, M. Soida and R.-J. Dettmar, Jagielonian University, Kracow, p. 231 

11. EnBHn, T., & Vogt, C. 2006, ApJ, in press 

12. Fabian, A.C. 1994, ARA& A, 32, 277 

13. Fabian, A.C, Mushotzky, R.F, Nulsen, P.E.J., & Peterson, J.R. 2001, MNRAS, 321, L20 

14. Fabian, A.C, Reynolds, C.S., Taylor, G.B. & Dunn, R.J. 2005, MNRAS, 363, 891 

15. Goldreich, P & Sridhar, S. 1995, ApJ, 438, 763 

16. Higdon J. C, 1984, ApJ, 285, 109 

17. Kota, J. & Jokipii, J. 2000, ApJ, 531, 1067 

18. Kulsrud R., 2004, Plasma Physics for Astrophysics, Princeton, NJ, Princeton University Press 

19. Lazarian, A. 2007, ApJ, 660, 173 

20. Lazai-ian, A. 2006, ApJ, 645, 25 

21. Lazai-ian, A. & Beresnyak, A. 2006, MNRAS, 373, 1195 

22. Lazai-ian, A., Vishniac, E. & Cho, J. 2004, ApJ, 603, 180 

23. Lazarian, A. & Vishniac, E. 1999, ApJ, 517, 700 

24. Malyshkin, L. & Kulsrud, R. 2001, ApJ, 549, 402 

25. Lesieur, M. 1990, Turbulence in fluids : stochastic and numerical modelling, 2nd. rev. ed. (Dordrecht; 
Kluwer) 

26. Narayan, R., & Medvedev M. 2001, ApJ, 562, LI 29 

27. Nusser, A., Silk, J. & Babul, A. 2006, astro-ph/0602566| 

28. Piffaretti, R. & Kaastra, J. 2006, astro-ph/0602376 

29. Rechester, A., & Rosenbluth, M. 1978, Phys. Rev. Lett., 40, 38 

30. Ruzkowski, M. & Begelman, M.C 2002, ApJ, 581, 223 

31. Schekochihin, A. & Cowley, S. 2006, asti-o-ph/0601246 

32. Shebalin J.V, Matthaeus W.H., Montgomery D.C, 1983, J. Plasma Phys., 29, 525 

33. Smith, R. & Cox, D. 2001, ApJS, 134, 283 

34. Voigt, L.M. & Fabian, A.C. 2004, MNRAS, 347, 1 130 

35. Webb, G., Zank, G., Kaghashvili, E., & Roux J. 2006, ApJ, 651, 211 



